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Abstract 

In this paper we use AdS/CFT ideas in conjunction with insights from finite 
temperature real-time field theory formalism to compute 3-point correlators of N=^ 
super Yang-Mills operators, in real time and at finite temperature. To this end, we 
propose that the gravity field action is integrated only over the right and left quad- 
rants of the Penrose diagram of the Anti de Sitter-Schwarzschild background, with a 
CN ' relative sign between the two terms. For concreteness we consider the case of a scalar 

Q"^ ' field in the black hole background. Using the scalar field Schwinger-Keldysh bulk-to- 



boundary propagators, we give the general expression of a 3-point real-time Green's 
correlator. We then note that this particular prescription amounts to adapting the 
finite-temperature analog of Veltman's circling rules to tree-level Witten diagrams. 



o . 

. and comment on the retarded and Feynman scalar bulk-to-boundary propagators. 

We subject our prescription to several checks: KMS identities, the largest time 
, ■ equation and the zero-temperature limit. When specializing to a particular retarded 

^ ■ (causal) 3-point function, we find a very simple answer: the momentum-space cor- 

^ relator is given by three causal (two advanced and one retarded) bulk-to-boundary 

propagators, meeting at a vertex point which is integrated from spatial infinity to 
the horizon only. This result is expected based on analyticity, since the retarded 
n-point functions are obtained by analytic continuation from the imaginary time 
Green's function, and based on causality considerations. 



^E-mail addresses: eb4df,dv3h,cw2an, parnold@virginia.edu 



1 Introduction and Summary 



In recent years, there has been a great deal of interest in applying AdS / CFT methods to 
the study of strongly-coupled quark-gluon plasmas created at the Relativistic Heavy Ion 
Collider (RHIC) at Brookhaven National Laboratory. This interest stems from the fact 
that some features of strongly-coupled quark-gluon plasmas are captured by an analyti- 
cally tractable example of a gauge theory plasma at strong coupling: A/'=4 SU{Nc) super 
Yang-Mills theory at finite temperature. In the deconfined phase, this theory is expected 
to qualitatively resemble high-temperature QCD (see e.g. [1] for a discussion of this). In 
the limit of large Nc and large 't Hooft coupling, A/'=4 super Yang-Mills is holographi- 
cally dual to weakly-coupled string theory on a black hole spacetime which asymptotes 
to the product of 5-dimensional Anti-de Sitter space (AdS) with a 5-dimensional sphere 
[2, 3, 4, 5]. Using the AdS/CFT correspondence, many dynamical properties of the A/'=4 
strongly coupled plasma have been investigated: e.g. transport coefficients such as vis- 
cosity were related to a computation of real-time 2-point correlators in the black hole 
background (see [6] and references within). The prescription for computing 2-point real- 
time correlators goes back to Son and Starinets [7] and to Son and Herzog [8] (for more 
recent work see [9]). On the other hand, very little is known about computing 3-point and 
higher correlators of operators in real time and at finite temperature. In fact, this paper 
gives the very first concrete expression for real-time finite-temperature 3-point correlators. 

According to AdS/CFT, one can compute a 3-point correlator in super Yang-Mills 
by computing the analog of a Feynman diagram (Witten diagram) in the Anti-de Sitter- 
Schwarzschild (AdS-S) background. This diagram includes three bulk-to-boundary prop- 
agators for the gravity field which couples to the super Yang-Mills operator under con- 
sideration, and these propagators meet at a vertex point in the bulk. The location of the 
vertex must be integrated over. For real-time (Minkowski signature) calculations, a possi- 
ble source of confusion concerns the range of integration. Should one integrate only up to 
the horizon in Schwarzschild coordinates? Or should one perhaps think about the Penrose 
diagram of the AdS-S space and integrate over the right and left quadrants only? One 
could instead integrate over the whole diagram (including the past and future quadrants), 
but then one would need to specify how to carry out the integration at the singularities. 
Or maybe one has to use an entirely different gravitational background by gluing different 
geometries in such a way that one constructs a gravitational analogue of the Schwinger- 
Keldysh contour. Various authors have proposed different ways of addressing this issue 
[8], [9, 10]. We will resolve this integration ambiguity by extending the work of Son and 
Herzog and using many insights from the formalism of real-time finite-temperature field 
theory, offering an alternative to other approaches in the literature which is quite simple 
to implement, and which yields consistent results. 

For simplicity, we consider operators which couple to (massless) scalar fields in the 
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AdS-S background, but we expect that our results generalize easily to other spins. We 
first construct the Schwinger-Keldysh scalar bulk-to-boundary propagator. This is a 2 x 2 
matrix whose indices correspond to the location of the boundary point (either the right 
(R) or left (L) boundary) and to the location of the bulk vertex point (either in the R or 
L quadrants). The minimally coupled scalar action is defined to be integrated over only 
the R and L quadrants, with a relative sign between the two bulk integrals. This is along 
the lines of Frolov and Martinez [11], who in turn were inspired by Israel's thermofield 
formalism [12]. It turns out that the relative sign between the two contributions precisely 
matches the relative sign between the contributions of the physical and doubler fields in 
real-time finite-temperature field theory. We note, following Kobes and Semenoff [13], that 
the Schwinger-Keldysh propagator can be re-interpreted in terms of circling rules (that 
is, the entries of the Schwinger-Keldysh propagator can be related to the Feynman and 
Wightman propagators, which are represented diagrammatically in terms of Veltman's 
circling rules [14, 15]). Thus the R-L prescription can be traded for circling rules in the 
AdS-S background, and with the integration over the bulk vertex only up to the horizon. 
The various real-time 3-point functions are then computed diagrammatically by placing 
circles around the boundary and bulk vertex points, and summing the appropriate number 
of diagrams. For example, a time-ordered product of three operators in real time and 
finite temperature is computed by adding two gravity diagrams: one with three boundary 
vertices of type 1 (uncircled) and three bulk-to-boundary propagators meeting at a bulk 
vertex which is either of type 1 (uncircled) or 2 (circled). As we have already mentioned, 
the bulk vertex is integrated only up to the horizon. 

We subject our prescription to several checks: Kubo- Martin- Schwinger identities, the 
largest time equation, and the zero-temperature limit. Moreover, by focusing on a special 
3-point real-time correlator which is retarded, we can make contact with the 3-point 
correlator obtained in imaginary time (Euclidean signature). The latter can be computed 
straightforwardly using Witten-type diagrams in Euclidean AdS-S, without encountering 
any of the subtleties we discussed before: Euclidean AdS-S has no singularities (it caps 
at the horizon) and has only a single boundary, just like Euclidean AdS. The retarded 
3-point function is related to the imaginary-time 3-point Green's function by analytic 
continuation in frequency. Our more general prescription for computing generic real- 
time finite-temperature 3-point correlators yields a retarded 3-point function which has 
this property. In the past, analyticity arguments were used in the context of AdS/CFT 
by Gubser et al. [16] and by Iqbal and Liu [17, 18] to obtain the real-time 2-point 
correlators. Lastly, using the causal nature of the retarded 3-point, we can understand 
why the integration region over the bulk vertex is only up to the horizon, that is, why 
one needs to integrate only over the R region of the Penrose diagram (similar causality 
arguments were used by Caron-Huot and Saremi [19] when computing one-loop gravity 
corrections to the retarded 2-point correlator). 
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The paper is organized as follows: In Section 2, we first derive the real-time 3-point 
correlators at zero temperature. This is done so that we have a reference point when asking 
the question: is the finite-temperature prescription compatible with the zero-temperature 
result? In our approach we start from the position-space field theory correlators and 
Fourier-transform. Amusingly, in this procedure, the AdS radial coordinate arises as 
a Schwinger parameter. Next we find that the real-time correlators are obtained by 
integrating over a region of the AdS space, the Poincare patch. This is in contrast to the 
more familiar AdS/CFT story in Euclidean signature, where the AdS tree-level diagrams 
are obtained by integrating over the whole AdS space. In particular, the final answer 
for the momentum-space retarded 3-point correlator is expressed as the product of three 
causal (two advanced and one retarded) bulk-to-boundary AdS propagators, integrated 
over the position of the bulk vertex only over the Poincare patch. 

Then, in Section 3 we review some of the basic elements of real-time finite-temperature 
field theory formalism. In particular we take note of an observation made by Kobes [20] 
that in real-time formalism, the retarded n-point function, which he proceeds to define, 
is the real-time correlator which is obtained by analytic continuation in frequency space 
from the imaginary-time finite-temperature correlator. 

In Section 4, we spend some time reviewing the construction of the bulk-to-boundary 
propagators. In particular, the retarded bulk-to-boundary scalar propagator is the one 
which behaves like an incoming wave at the horizon [19]. We give a self-contained expo- 
sition of its expression in terms of Heun's functions, and explore its analytic properties, 
connection with the Euclidean signature bulk-to-boundary propagator, causal properties 
and zero-temperature limit. Then we notice that the Feynman propagator in curved space 
is associated with a choice of vacuum. We construct the Feynman Green's function from 
the retarded Green's function by borrowing the definition of a thermal Feynman Green's 
function from finite-temperature field theory: 

GpiE, P, u) = ReGniE, P, u) + zIm^R(E, P, n)coth(/3E/2), (1.1) 

where u is related to the radial/holographic direction, and f3 is the inverse tempera- 
ture. We then note that the Schwinger- Keldysh bulk-to-boundary propagator can be 
re-expressed in terms of circling rules propagators [14, 15] in a way entirely analogous to 
real-time finite-temperature field theory [13]. 

We begin Section 5 by reviewing the real-time finite-temperature 2-point correlators 
from AdS/CFT. Our definition for the scalar field action is in accord with the R-L pre- 
scription. This choice is such that by evaluating the on-shell AdS-S scalar kinetic action 
and picking up the R and L boundary terms, one gets the 2-point Schwinger-Keldysh field 
theory correlator. Next, we use the same R-L prescription to evaluate the 3-point real- 
time finite-temperature correlators. We have already advertised the result: it is consistent 
with the known identities obeyed by the real-time finite-temperature correlators, and with 
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the zero-temperature limit. Then, by focusing on the retarded finite-temperature 3-point 
correlator, we find that although numerous terms contribute, the answer is quite simple: 
it can be viewed as arising from a single Witten-type diagram containing three causal 
propagators (two advanced and one retarded) joined at a bulk vertex which is integrated 
only up to the horizon. This is what one expects to get based on considerations of the 
zero-temperature limit and analytic continuation from Euclidean signature. Moreover, for 
a causal 3-point function, one expects that one has to integrate only over a causal bulk 
region determined by the boundary points. Since our result gives the momentum-space 
3-point retarded correlator, the bulk integration can only be over the R quadrant, which 
is a maximal causal diamond. We view all of these results as consistency checks of our 
prescription for the computation of the real-time finite-temperature 3-point functions, 
which can also be summed up by the use of circling rules, and a bulk integration region 
from spatial infinity up to the black hole horizon. 

The technical details are relegated to Appendices. In Appendix A, for completeness, 
we give the momentum-space zero-temperature 2-point CFT correlators. In Appendix 
B, we obtain the real-time momentum-space zero-temperature 3-point functions in terms 
of Witten-type diagrams by employing a reverse-engineering procedure, namely starting 
from the CFT correlators and re-packaging them as Witten diagrams. In Appendix C, we 
discuss the zero-temperature retarded bulk-to-boundary scalar propagators. In Appendix 
D, we give the closed-form expression of the zero-temperature retarded 3-point function 
for CFT operators with conformal dimension A = 2. Lastly, in Appendix E we give some 
background material on the Heun function and its various local representations. 

2 From field theory T=0 real-time correlators to AdS 
Witten diagrams 

2.1 Real-Time CFT correlators 

To set our notation in a simplified setting, let us begin with a massless scalar field theory. 
In Minkowski signature (that is, real time), one encounters different Green's functions. 
Besides the Feynman propagator. 
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Figure 1: Circling rules propagators 



where are Wightman's functions 



D^ix) = (O|0(O)0(x)|O) 
1 1 



(27r)2 [-(t + ie)2+f2] 

D+{x) = = {0\(j){x)(j){0)\0) 
1 1 



(27r)2 [-(t-ze)2 + f2]' 
one defines the causal retarded and advanced propagators, 

Dnix) = ^(t)([0(x),0(O)]) 

= Df{x)-D~{x) 
i 



da{x) ^ -e{-t)mx),m]) 

= Df{x)-D+{x) 



(2.2) 
(2.3) 



-^(t)5(-t2 + f^) (2.4) 



-^(-t)5(-t^+x2). (2.5) 

The retarded/advanced propagators (and the higher n-point functions) can be obtained 
following Veltman's circling rules [14, 15]. According to these rules, each vertex can be 
circled or uncircled. The circling of a vertex brings a minus sign. The propagator between 
uncircled vertices is the usual Feynman propagator, Dp{x). The propagator between 
two circled vertices is the complex conjugate of the Feynman propagator, Dp{x). The 
propagator between a circled vertex and an uncircled one is D~{x), and the propagator 
between an uncircled vertex and a circled one is D~^{x) (see Fig.l). As a consequence 
of these circling rules, Veltman was able to formulate the largest time equation [14]. 
This is an algebraic identity stating that the sum of all diagrams obtained from a single 
Feynman diagram by placing circles around all vertices in all possible combinations, for 
a total of 2^ diagrams if there are n vertices, is zero. For example, one can easily verify 
that Dp- D- - D+ + D*p = 0. 
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The causal retarded n-point function where the largest time is associated with a certain 
vertex is computed by adding all diagrams, where each vertex can be circled or uncircled, 
with the exception of the vertex associated with the largest time which remains uncircled. 
For the 2-point, this indeed reduces to Dr{x) = Df{x) — D~{x). 

Our notation is such that the propagator will be consistently denoted by fonts of the 
letter D, while the n-point Green's functions will be denoted by fonts of the letter G. For 
example, the 2- and n-point Feynman Green's functions are given by 

iGf{x, x') = (r0(x)0(x')), i^-'Gf{xi,X2 . . . x„) = (r0(xi)0(x2) . . . (2.6) 

In a conformal field theory, the 2-point function is fixed by symmetry. Considering 
for simplicity scalar operators of the same conformal dimension A, the various real-time 
2-point functions are^: 

Gf{x) = -^{0\TOix)OmO) = -j^J ,2 ^ ■ ) ' (2-7) 



A 



(27r)2 \-(t-iey + x^ 

(2.8) 

A 



G+{x) = -«(0|O(x)O(0)|0) = 9{t)GFix) - ei-t)G*pix) 



G~{x) = -z{0\O{0)O{x)\0) = e{-t)GF{x)~e{t)GUx)= 

(2.9) 

Gr{x) = Gf{x) -G-{x) = e{t)(G+{x) -G-{x)\ (2.10) 



GAix) = GFix)-G+{x) = e{-t)\^G~{x)-G+{x)j. (2.11) 

The CFT Green's functions have the same ie prescription as the scalar Green's functions. 
For completeness, we give the momentum-space 2-point functions in Appendix A. An 
important observation which carries through, as will see when discussing the retarded 
momentum-space propagator of scalar fields in a black hole background, is that the an- 
alytic continuation of the Euclidean signature 2-point function, with E — —i{E± ie) 
yields the retarded/advanced 2-point correlators [16] [17, 18]. 

Moving on to 3-point functions, these are also fixed by conformal symmetry (up to an 



^The A/'=4 super Yang-Mills theory 2-point functions have an additional overall constant 8^^^ r(A-2) 
which was stripped off from the subsequent formulae, together with a factor which is dependent on the 
gauge group and representation. 
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overall constant, which was set to 1). For example 

G^(Xi,X2,X3) = (-2)2(0|rO(Xi)O(x2)O(x3)|0) 



1 



_f2 , ^ 
''12 ' -^12 



% + ^23 + + zel 



A/2 

(2.12) 



The non-time-ordered product of three scalar operators is 



H2 + ^12 



ieti2) (-4 + ^23 



^£^23) (-^ii + ^ii - i(^hi) 



A/2 



(2.13) 



where the ie prescription follows the same rules as for the 2-point functions, namely the 
time coordinate of an operator insertion is greater by ie than the time coordinate of any 
other operator insertion to the right of it [21] (for a recent AdS/CFT paper making use 
of this ie prescription see [22]). 

Knowing these correlators means that we can obtain the retarded 3-point function for 
which X3 has the largest time by using the analytic continuations which follow from the 
circling rules (alternatively, see [23]): 



Gr{Xi, X2; X3) — 9{t^i)d{ti2) I G312 — (-^132 + 6*213 ~ 6*231 



= {-iYe{h{)e{t^2 

+ (-^)'^(^32)^(^2l 



+^(^32)^(^21) ( 6*321 — 6*231 + 6*123 ~ 6*132 
1 1 



{xl^ - letuY (xi3 + ^6^23)"^ 

1 1 
{XI2 + ^6^12)^ (xii - ietsiY 



c.c. 



c.c. 



(xii - letsiY (x|i + letsiY 



[^3 + ^(^hsY 



(xl^ - iet23Y 



(2.14) 



where we defined for brevity 
and Gabc = {—iY{0{xA)0{xB)0{xc)) are the non-time-ordered correlators. 



2 ' 



(2.15) 



2.2 Lessons for real-time AdS/CFT from CFT correlators 

Before we discuss the real-time correlators, we make a detour in Euclidean space, and 
show how one could discover Witten's AdS/CFT diagrams starting from the 3-point 
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momentum-space CFT correlator. We begin by Fourier-transforming the position-space 
3-point correlator. To perform the integrals, we introduce three Schwinger parameters. 
The integrals over the Euclidean-signature coordinates x'^ are of Guassian type, and easily 
performed. To deal with the intermediate result, we need to introduce one more Schwinger 
parameter, z. As shown in detail in Appendix B, the final answer is: 



i'^GEih, k2, k 



3 J 



4 2 ^ 
ol0-3A-3e^^,^2^4-26 ^ ~ ^(^1 + ^2 + ^3) dz /7^. A-^:^ tv- / /T2 



(2.16) 



The expression in (2.16) is UV-divergent, and to regularize it we employed here dimen- 
sional regularization. It is amusing to notice that the momentum-space 3-point CFT 
correlator has been reassembled as an integral over the AdS radial coordinate z of the 
product of three bulk-to-boundary AdS scalar propagators. This is precisely a 3-point 
function Witten diagram [3]. With a bit of hindsight, from the measure factor and from 
the form of the propagators in (2.16), the AdS metric can be reconstructed to yield 



dsl = + (2.17) 



In Euclidean signature the Poincare coordinates of (2.17) cover the whole AdS space. 

We should also point out that as a result of our use of dimensional regularization in 
computing (2.16), the AdS space has dimension d + 1, with d = 4 — 2e. The scalar field 
with mass m and d-dimensional momentum /c^ in {d + l)-dimensional Euclidean AdS 
space has a bulk-to-boundary propagator 

{^)^-iziK^_d{^z), (2.18) 

where A is the conformal dimension of the CFT operator which couples to mass m scalar 
field: 

A = ]^{d + y/d?T^). (2.19) 

A different regularization which is more commonly used in AdS / CFT calculations im- 
poses boundary conditions aX z = zb ^ I [24]. Then, the bulk-to-boundary propagators 
which have a 5-function support at in position space, have the following expression in 
momentum space: 

z^K^_d{^z) 

. (2.20) 



2 



zhKA-^WkfzB) 
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Next, we would like to compute the Fourier transform of (2.14) to get the momentum 
space retarded 3-point correlator. The reason to perform this calculation is to find out 
how the Witten diagrams look like in real-time. In particular we want to find out which 
is the region of the Minkowski-signature AdS space where the tree-level gravity diagrams 
need to be integrated over. This becomes relevant later, in Section 5, when we will 
inquire whether the real-time finite-temperature 3-point correlators have the correct zero- 
temperature limit. 

As we will see, the answer is quite natural: since we want to identify the 4-dimensional 
coordinates in the field theory and in the holographic dual, we will be using the Poincare 
parametrization of the AdS metric (2.17) with Minkowski signature. This also comes out 
naturally from the reverse engineering perspective that we are following: as in (2.16), this 
identification is built-in. It follows that we will not be integrating over the whole AdS 
space, since (2.17) now covers just half of AdS. As for the question of whether in real-time 
AdS/CFT we will be integrating over the whole range of the radial coordinate z, we can 
anticipate that the answer will be affirmative, if the radial coordinate enters, as it did 
before, as a Schwinger parameter. Lastly, we might expect that the retarded 3-point is 
obtained by performing an analytic continuation uj —i{u! + ie) of the Euclidean result, 
and that perhaps this analytic continuation is allowed under the integral (2.16). The 
answer to this second question is again affirmative. 

The most straightforward proof of our previous statements requires that we choose a 
certain momentum kinematics for the momentum-space correlator: if the time associated 
with the spacetime point xs is the largest of xi, X2, X3, we will take the momenta Pi 2 3 to 
be such that £'1,2 < 0, -Es > 0. With these kinematics, we can show that our retarded 
3-point correlator in momentum space is given by the Fourier transform of Gf{xi, X2, X3). 
To see this we only need to repackage the retarded 3-point in position space as 



i^GR{Xi,X2]X3) 





1 


1 


1 


''12 


h xf 2 


+ ie —^23 "1" -^23 "1" ~^31 


+ xl^+ie) 




1 


1 


1 


-(tl2 


-^e) 


2 + X12 ~^23 "1" -^3 + ~ 


(t3i + leY + xli 




1 


1 


1 


-(tl2 


+ ie) 


^ +^?2 -{hs - ie)^ + 


— t|i + xli + ie 




1 


1 


1 


''12 


h X12 


- -(% - + 3^23 ~ 


{hi + ief + 



A/2 



A/2 



A/2 

(2.21) 

and then write the Fourier transform in terms of an integral over a "loop momentum" and 
use the results given in Appendix A for the Fourier transform of each of the three factors 
composing each one of the four terms in (2.21). The energy step-functions that enter the 
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last three terms are incompatible with the kinematics that we have chosen. That is, with 
these kinematics, the first term alone accounts for the retarded 3-point function. 

As shown in Appendix B, the final answer for the momentum-space retarded 3-point 
correlator given in (2.21) is: 

Gn{Pi,P2;Ps) = 23A-6r(A/2)3r(3A/2-2) ^^^^+^^+^^Vo ^ 

z\-{E, - zeY+f,)^-'K^.,{z^-{E2 - lef + p^) 
z'{-{E, + ie)f + f,)^-'K^.,{z^-{E, + lef + pf) . 

(2.22) 

Each of the three factors present in the previous formula (e.g. z^{—{E-i+ieY+p^)^^^ 
KA-2{z^y—{E3 + ieY + p|)) is a causal (retarded in this case) bulk-to-boundary scalar 
propagator in the AdS background (see Appendix C). The retarded bulk-to-boundary 
propagator in AdS can be obtained by analytic continuation, E ^ —i{E + ie), from the 
Euclidean propagator (2.18). 

The integral on the right-hand side of (2.22) is not convergent for A > 2 (the case A = 
2 yields a convergent integral whose closed-form analytic expression is given in Appendix 
D). This is dealt with in the usual manner: the integrals could have been regularized 
using dimensional regularization, as we did before in manipulating the Euclidean 3-point 
correlator, or one employs boundary cut-off regularization. In this case, the integral over 
the radial coordinate z is cut-off at zb <C 1, and the propagators are replaced by functions 
which are normalized at the boundary as in (2.20). 



3 Field theory T^O real-time formalism 

In this section we review a few of the fundamental notions and definitions of the field 
theory real-time finite-temperature formalism. 

In real-time formalism one distinguishes between physical fields 

0i(x) = 0(f,t) (3.1) 

and doubler fields 

02 = 0(x, t - id) (3.2) 
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where a is the arbitrary parameter of the Schwinger-Keldysh contour. Correspondingly, 
the real-time n-point Green's functions are defined as 



Ga,a,...aS^.2,...n) ^ (-zf-i (Tp^,, (1)0,, (2) . . . (n));3, (3.3) 

where ai, . . . a„ = 1, 2. By 7p we denoted the time ordering along the Schwinger-Keldysh 
contour, where the 0i(x) fields are time-ordered, the 02 (a^) doubler fields are anti-time 
ordered, and lastly with any field 02 being defined to have a larger path-time than any 
other field 0i. In (3.3) the brackets (...)/? denote the thermal average over all Hamilto- 
nian eigenstates, each being weighted by the Maxwell-Boltzmann factor. The generating 
functional for the Green's functions is 



Z[ Ji, .h] = {Tp exp (^z J d^x{Jih - J202) J (3.4) 

and so 

G.,.,....„(l, 2, . . . n) = .(-1)^--"' . nuTlt/'U ( Y (3.5) 
If CT = /3/2, then one can prove that 

^*aia2...a„ = (~1)" ^Gaia2...a„ (3.6) 

where 1 = 2 and 2 = 1. These are the Kubo-Martin-Schwinger identities. 

Consider now a massless scalar field theory at finite temperature in the real-time 
formalism. The propagator is a 2 x 2 matrix 

Duip) DM\ ^[ + n{\E\)27rS{p') {e{-E) + n{\E\))27^5{p^)eP^I^' 

DM DM) \{0{E) + n(|E|))27r5(/)e-/5^/2 + n(|£;|)27r5(p2) 

(3.7) 

where 

"(1^1) = (^-^^ 

and /3 = 1/T is the inverse temperature. The off-diagonal components of the Schwinger- 
Keldysh propagator (3.7) are related to the finite-temperature Wightman functions D^{p): 

Di2 = exp{f3E/2)D-, D^i = exp{-PE/2)D+ , (3.9) 

while the diagonal components are related to the finite-temperature Feynman propagator: 

D^^ = Dp, D22 = D*p. (3.10) 
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In more generality, if the parameter a of the Schwinger-Keldysh contour is not equal to 
/3/2 as in the previous formulae, then 

Du = Df, D22 = D*f, Du{p)=exp{aE)D~{p), D^iip) = exp{-aE)D+{p). 

(3.11) 

With the choice a = (3/2, the Schwinger-Keldysh propagator is a symmetric matrix, as 
in (3.7). As it turns out, the finite-temperature computations performed via AdS/CFT 
following [8] yield 2-point correlators with a = /3/2. 

The matrix propagator (3.7) contains the propagators between all types of vertices: 
physical field vertices being labeled 1, and doubler field vertices, labeled by 2. A real-time 
Feynman Gn...! n-point function is computed diagrammatically by summing all distinct 
diagrams obtained by letting the vertices not connected to external lines be of either type, 
1 or 2. The vertices connected to external lines are taken to be of type 1. This turns out 
to be equivalent to just using the circling rules, and the Dp, Dp, D^ propagators between 
circled/uncircled vertices, with only the internal vertices being of either circled/uncircled 
type [13]. Moreover, the largest time equation of Veltman still holds at finite temperature. 

These observations led Kobes [20] to define a "causal" n-point function, by adding all 
diagrams with a vertex circled or uncircled, with the exception of the vertex associated 
with "the largest time" , which remains uncircled. This is very much the same prescription 
used at zero temperature. It is this retarded n-point function which is obtained by the 
analytic continuation of the n-point function computed in imaginary-time formalism [20]^. 

For example, a causal/retarded 3-point function, with the outgoing momentum r = 
— {p + q) vertex corresponding to "the largest time" is equal to 

Gn{q,p; r) = Gm - e-^^^/^G^^i - e-^^'/^G^n + e^^^/^G^^i- (3.12) 
In addition to the identities which follow from (3.6) 

111 — <-^222' '-^121 — *-^212' '^112 — <-^221' *-^211 — '^■122; (^O.iOj 

the largest time equation yields one more identity 

Q = G,,,-e-P''^''G,2i-e-f'''^''G2ii+e^''^''G22i+e^''^''G2i^^^^ 

(3.14) 

which, using (3.13), can be re-written as 

— sinh(a;j,7r)Re(j>'ii2(q',p) = smh.{(jjq'K)YieG2ii{q,p) + sinh(ci;p7r)ReGi2i(q',p), 
-ImGiii(g,p) + cosh(ci;p7r)ImGi2i(g,p) + cosh(cUq7r)ImG2ii(g,p) + cosh(u;,.7r)ImGii2(g,p) 

(3 

■^The other real-time finite-temperature correlators do not enjoy such a simple relationship with the 
imaginary-time correlator. 
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Substituting (3.15) into (3.12) leads to a simpler expression for the causal 3-point function 
[20]: 

T. ^ / ^ T. ^ smh{(3Eg/2)^ ^ sinh(/3Ep/2)^ ^ 

sinn(p_fcr/2) sinn(p_fcf,/2) 

lmGn{q,p;r) = - tanh(/3E./2) f ImGni + '^^i^^MG,,, + '^l[ifih ^G,n 

(3.17) 



4 Real-time AdS-S propagators and Circling rules 

Some of the results presented in this section are known in the hterature [19]. However 
we cound not find references giving a comprehensive picture of the material included 
here. Also since the literature is not conclusive in terms of how to best formulate real- 
time AdS/CFT [25, 26, 27, 7, 8, 9], we have decided to spend some time discussing 
the various bulk-to-boundary propagators (retarded, Feynman, etc..) in AdS-S. These 
propagators will consitute the basis of a diagrammatic expansion in terms of tree-level 
gravity diagrams, modeled on the finite-temperature real-time field theory formalism. 
The AdS-S retarded propagator is the one which behaves like an incoming wave at the 
horizon. We discuss its analytic properties, causal nature and zero-temperature limit. 
The Feynman propagator which corresponds to the thermal "Kruskal" vacuum is obtained 
from the retarded propagator using (1.1), which is typical of finite-temperature systems. 
Next we formulate the AdS-S circling rules as the curved space counterpart of the same 
rules introduced in the previous sections. This will help us develop a better perspective 
on the prescription for computing the real-time finite-temperature correlators in Section 
5. 



4.1 The AdS-Schwarzschild Geometry 

The Anti de Sitter-Schwarzschild (AdS-S) black hole (times a five-dimenional sphere 5"^) 
is the holographic dual of the finite-temperature A/'=4 super Yang-Mills theory with gauge 
group SU{Nc), in the limit iVc ^ 1, and in the deconfined (high-temperature) phase [5]. 
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The AdS-S metric is usually written as 



fir) 



1 - 



dr"^ + R^dVt. 



dnl ,, 



r 



(4.1) 



(4.2) 



where the position of the black hole horizon is at tq and the asymptotic region of the 
black hole geometry is at r = oo (or z = 0). The black hole Hawking temperature is 



(4.3) 



A better suited choice of coordinates is 



u 



u 



ds^ - ^ 



1 — , , dx^ 

dt^ + 

u u 



2nn2 p2 

^ ^H^.^_^l_u')dt' + dx') + 



4(1 - ^2)^2 

R' 



du 



:du . 



(4.4) 



u ' ' ' ' 4(1 — ^2)^2 

In these new coordinates, the horizon is at m = 1, the singularity is at u = oo and the 
asymptotic region is at u = 0. 

It has been conjectured that finite-temperature correlators should be related via AdS / CFT 
to n-point functions of supergravity living in the maximally extended AdS-S geometry 
[8]. Therefore, we will begin by switching to global coordinates which describe the ex- 
tended geometry. We will look for coordinates tK, xk which are the analog of Kruskal 
coordinates for Schwarzschild black holes in flat space. These can be found by starting 
with the ansatz: 



tx = TZiu) sinh {27iTh t) , xk = Tiiu) cosh {2tiTh t) 
In terms of tx and xk-, we demand that the metric have the form 



ds^ = W{u) [-dt\ + dx\] + 



i^TuRf,,, 



dx' + R'dn 



u 



5- 



Matching this metric to our black hole metric requires 

1 



7^(n) 
W{u) 



exp 



tan 



log 



1 - 



f{n) 
4u7^2(n)■ 



(4.5) 
(4.6) 

(4.7) 
(4.8) 
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Figure 2: The AdS-S Penrose diagram 



Since xj^ ~ tj^ = '^(^)^; the event horizon is described by two intersecting hnes in the 
tx/xK plane: 

event horizon: x\ - t\ = n^{l) = (4.9) 

In the above expression for TZ{u), we have also taken the liberty of setting to zero an 
integration constant which would appear as an overall rescaling of TZ. This fixes the 
distance between the origin [tx, xk) = (0, 0) and the AdS boundary to be 1: 

AdS boundary: - 4 = ^(0)^ = 1- (4.10) 

This choice also fixes the distance between the origin and the black hole singularity: 

singularity: x| - = 7^^(oo) = -e"''. (4.11) 

It is apparent from the above expressions for the event horizon, AdS boundary and black 
hole singularity that, in terms of and xk, there are two event horizons, two AdS 
boundaries, and two singularities, if xk is analytically continued to negative values. There 
is no obstacle to performing this analytic continuation, and the result is the extended AdS- 
S geometry. In the extended geometry, tx and xk have the ranges tj^ G (—00,00) and 
Xk G (—00, 00). The Penrose diagram of the AdS-S geometry is given in Fig. 2^. The 
fact that there are two AdS boundaries is particularly important for the computation of 
time-dependent correlators. 

In what follows we will ignore the fluctuations in the compact directions, and 
consider only gravity fluctuations in the 5-dimensional AdS-Scharzschild geometry (4.4). 

^We took the artistic license to draw a square Penrose diagram, but this is accurate only for the 
3-dimensional AdS-S black hole. 
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4.2 The retarded bulk-to-boundary scalar propagator in AdS-S 
4.2.1 Minkowski signature 

Consider a minimally coupled massless scalar propagating in the black hole background. 
To terms quadratic in fields, its equation of motion is 

d,{^gg^'d,cl>{x^,u)) = Q. (4.12) 

After Fourier-transforming along the coordinates to momentum space, the scalar field 
equation of motion becomes 



P"-4^,F^U—^^-^f^^F = ,, (4.13) 

M(1 — M^j \U[\ — W^Y M(1 — M 

where (j){p^,u) = F{u,p,u)(j)o{p'^), 0o(p^) is the boundary value of the scalar field, and 
prime denotes differentiation with respect to the bulk coordinate u. The variables a;,p 
are dimensionless quantities defined as 

E ^ P , , 

^ = — 7^,P= — 7^, (4.14) 
2nTH 2txTh 

where E, P are the energy and momentum associated with the scalar field modes. The 
norm of the dimensionless spatial momentum is denoted by The bulk-to-boundary 
propagator is further normalized to 1 at the boundary: F{u = 0) = 1. The 5-dimensional 
AdS-S bulk-to-boundary propagators admit analytic expressions in terms of Heun's func- 
tions (see also [7], [28]). 

We begin by making the substitution 

F = {l-u)^{l + u)^H. (4.15) 

Then, H obeys the differential equation 

H"+(-l + Zi±i^ + 1±^)h'+ ^<"' ^ IP^'' 7 + + '"'"g = 0. (4.16) 

\ u 1 — u l+u J 2u{l — w^) 

which is of Heun type (see Appendix E), with parameters^ 

d=-l, q = u^- \p\^ - ^^w, a = P = ^—^u, 7 = -1, 5 = l-iuj. (4.17) 



function. In this case, one starts with the ansatz F = u^{l — u)~2-(l + u)t: where \ = ^ oy 2 — — 
The Heun parameters wih then also depend on the conformal dimension A. 



The bulk-to-boundary propagator for a massive scalar field can also be expressed in terms of a Heun's 

|or2-#. 
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The two independent solutions of (4.16) are 



Hi(u) = lim 

^ 7-5>-l 



Hl{-l,q,a,f3,^,S;u), 



1 _ \^2^2 



-u 



1-7 



Hl{-l,q\a\P\i,5-u) 



2 7 + 1 

q' = uj'- - - (7 - 1)(^(1 + 



(4.18) 



1 — i 

a' = (3' = 00 - 7 + 1,7' = 2- 7 

and, 

H2iu) = u^Hli-1, q', a', (3', 7', S; u), (4.19) 

where in the definition of Hi 7 is kept arbitrary prior to taking the hmit (see equation 
(E.5) and related discussion in Appendix E). In H2, 7 can be set to —1 directly. For 
future reference we give the small u expansions of Hi{u) and H2{u): 

Hi{u) = 1 + ciu + {c2\n{u) + C2)u^ + . . . 

H2{u) = u^ + ... (4.20) 



where 



Ci=q = —u + u^ - IpI^ 

ju;' - \p\r 

C2 = |(l-9-c^(l + 0). (4.21) 

Next, we focus on constructing the solution of (4.16) which corresponds to a purely 
incoming wave at the horizon. This is the retarded bulk-to-boundary propagator in the 
black hole background [19]: 

F{u = 1) = incoming wave and F{u = 0) = 1^. (4.22) 



^For the case under consideration, that of a minimally coupled massless scalar, the bulk-to-boundary 
propagator remains finite at it = and therefore can be normalized at the boundary. However, the 3-point 
correlator is expressed as an integral over the radial coordinate u and has a divergent integrand at u = 
because of the measure factor. The integral is regularized usually by boundary cut-off rcgularization, in 
which case the bulk-to-boundary propagators are normalized at ub ^ 1 • 
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An incoming wave at the horizon (u = 1) behaves as (1 — u) , whereas the outgoing 
wave is (1 — u)^ . 

To this end, we make a change of variable, 

w = l-u, H{u)=H{w). (4.23) 

After substituting into (4.16), we arrive at another Heun equation for H. The incoming 
wave at the horizon corresponds to the solution 

H,{u) ^ m{2, -'-) + ^^co + \p\', ^uj, ^co, 1 - zco, -1; 1-u). (4.24) 

Finally, the bulk-to-boundary propagator corresponding to an incoming wave can be ex- 
pressed as 

F = B{l-uy'^{l+u)^H3{u), (4.25) 
where 5 is a normalization coefficient 

Alternatively, we can write 

F = (1 - uy'^il + u)^ (^Hi{u) + AH2{u)^ , (4.27) 
where the coefficient A is^ 



A can be solved for in the limit of small u and directly from (4.16), by writing 
H{u) = 1 — 'Hh{u) = 1 — {^uj + |p|^)/2 h{u) and solving for h{u) such that to order Ti 
the solution one finds is regular at the horizon and it is normalized to 1 at the boundary. 
This allows computing the terms linear and quadratic in m, to order Ti. Recalling that 
the terms quadratic in m in if have a coefficient C2 + A, we get 

A = ^^co + + 0(a;^ a;bf , bt). (4.29) 

Earlier we have advertised that F is the retarded bulk-to-boundary propagator. To 
better understand that is so, we will show that F is obtained by analytic continuation of 



^In practice, in numerical computations we have matched (4.25) and (4.27) at some value Um = 0.5 
where both (4.25) and (4.27) are within the radius of convergence of the corresponding Heun's functions. 
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the Euclidean signature propagator which is regular at the origin (i.e. at n = 1) for values 
of the Euchdean frequency such that Re(ci;^) > 0. Then we will argue that F is analytic 
in the upper half w-plane, and its Fourier transform to position space is a causal function, 
with support inside the future light cone. Lastly, we will show that the zero-temperature 
limit of the AdS-S retarded propagator is, as expected, another causal propagator, the 
AdS retarded bulk-to-boundary propagator. 



4.2.2 Euclidean signature 

Here we address the same problem of a bulk-to-boundary propagator, but in the Euclidean 
version of the AdS-S black hole: 

dsl = ^JIIe^ [(1 - u')dtl + dx'] + ^ du\ (4.30) 

The metric in (4.30) has a conical singularity at u = 1, unless the Euclidean time is 
periodic with period 1/Th- Then, the origin is a regular point, which corresponds to the 
black hole horizon in Minkowski signature. Also in the Euclideanized geometry, there is 
only a single (Euclidean) AdS boundary. One way to see this is to return to our earlier 
analysis of Kruskal coordinates. We could imagine writing a Euclidean analog of equation 
(4.5): 

tir = 7^(M)sin(27^T/^^s), xk = TZiu) cos {27iTh ts) . (4.31) 

In terms of these coordinates, we no longer get two disconnected boundaries and two 
horizons. Instead, the horizon is just a point [tx = xk = 0), and there is a single 
connected boundary which is just a circle of radius one: t"]^ + xj^ = 1. The black hole 
singularity is not part of the Euclidean geometry. 

The wave equation for a massless scalar in this geometry is 



n2 ,, . ,2 



n(l — n^) u{l — u 



2\2 



(w|; + (l-n')|A;n0 = O. (4.32) 



We are interested in positive frequency solutions that are regular at the origin {u = 1), 
so we set 

Fe{u) = 4^ = (1 - ur^/\l + uy^^/^Hiu, k) (4.33) 
0o(fc) 

where the dimensionless Euclidean 4- vector k'^ is k'^ = {ue, k). 0o(^) is the value of the 
scalar field on the boundary of Euclidean AdS-S. Given the periodicity of t^, we must 
conclude that ue is integer-valued. On the other hand, for the purpose of performing 
analytic continuation, we will allow ue to become complex- valued. 
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Plugging this ansatz into the equation for 0, we find that H satisfies the Heun equation: 

+ ' g^o. (4.34) 

\u u — 1 u — dj u[u — l)[u — a) 

with the various parameters given by 

l + i 

7 = -l, d=l + ujE, e = l + tujE, a = l3 = ^—UE, 

q = ^UJE-col-\k\\ d = -l. (4.35) 

It is easy to check that the condition e = a + /3 — 7 — (5 + lis satisfied. As in the previous 
section, we can write the solution to the Heun equation either in terms of Heun's functions 
defined in a neighborhood of the boundary m = or in terms of functions defined in a 
neighborhood of the horizon -u = 1. In the former case, we have 

FE = il- ur^/\l + + AH2{u)), (4.36) 

where we normalized the Euclidean propagator u = 0. Hi and H2 are two independent 
solutions to the Heun equation which are defined as in (4.18) and (4.19), and whose 
parameters are given in (4.35). 

Alternatively, we may express Fe in terms of a Heun's function defined near u = 1*: 

^ _ (1 - u)'^^/'^{l + uY'^^/'^Hlil -d,a(3- q; a, /3, 5, 7; 1 - u) 

HI{1 — d, af3 — q; a, /3, 6, 7; 1) 

By matching these two expressions for Fe, as in the previous section, we have 



(4.37) 



A 



(4.38) 



Using (4.37) we have generated plots for a wide range of the parameters ue, k"^ and u. 
We are only looking at the right half ue complex plane since the solution (4.37) gives a 



^Thc solution Fe which is regular at the origin for negative frequencies is 

(1 - + uY'^^^^Hljl -d,-q+{5- l)7d + (a - ^ + - ^ + 1); - (5 + 1, a - ^ + 1, 2 - J, 7; 1 ~ u) 

Hl{l -d,-q+{5- l)-fd + (a - (5 +!)(/?- (5 + 1); - (5 + 1, a - (5 + 1, 2 - (5, 7; 1) ' 

The analytic continuation of this Euclidean solution is the advanced propagator, which in momentum 
space is the complex conjugate of the retarded propagator F. 
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regular Euclidean propagator only for Iie{uE) > 0. From these plots (a couple of which 
are being shown in Figure 3) we have determined that is a smooth function with no 
infinite singularities for lie (ue) > 0. This is the expected behavior of an analytic function 
of cj£; in this region of the complex plane. The retarded propagator which is obtained 



Im(F£) 





Im((j£) 



Re(6L»£) 



-10 -10 

Figure 3: Re and Im parts of the Euclidean bulk-to-boundary propagator Fe, with 

u = 0.5 and k"^ = 1. 



from the Euclidean propagator by the usual analytic continuation, 

F{uj,p,u) = FE{-i{uj + ie)),p,u), 



(4.39) 



enjoys the same properties (smooth, free of infinite singularities) in the upper half Minkowski 
frequency plane Im(cj) > 0. 



4.2.3 Causality 

We will next argue that the Fourier transform of the finite-temperature retarded propa- 
gator, F{uj,p,u), is causal, i.e. has support inside the forward light cone. This is closely 
related to the analytic properties of F. Our analysis will be restricted for simplicity to the 
case when two spacetime points are separated in t and u but have the same coordinate 
f = 0. 

First, we construct the null AdS-S geodesies at x = 0: 







1 — rf 



u 



i?2 



4(1 - m2)m2 



(4.40) 
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with solution 



tiu) = nn 



If / 1, 

— arctan(V?i) + 2 log ( Y" 



t{u). (4.41) 



2nTH 

Values of t lying outside the forward light cone obey the condition t < t{u). 

Next we infer the behaviour of F{uj,p,u) at large frequency {u ^ 1) and fixed u^: 

F = y^{Lo + ze)3/2u3/4e-^e'"^(") + . . . , (4.42) 



Plots of the exact F along with this leading order behavior are shown in Fig. 4. When 
Fourier-transforming to position space, the leading term in (4.42) yields 

3y2^M3/4 I 

-9{t-t{u))- — (4.43) 



8 ' ' " {2^lTHt-T{u)fl^'' 
wii\l support inside the forward light cone. The next subleading term in (4.42) behaves 
like {(jj + ieY/'^e'^^'^^'^\ Its Fourier transform also has support inside the forward light cone. 
In order to analyze the behavior of the remaining contributions to F, we consider the 
following integral: 



oo 



— / due~'^' 
2vry_oo 

(4.44) 

where f{u,p, u) refers to the first sub-leading term in (4.42). For those values of t which 
lie outside the light cone {t{u) — t > 0), we may close the contour in the upper half u 
plane. The integrand has no poles or cuts in this region, and so the integral vanishes. 
Since the retarded propagator is causal for a; = 0, and since the zero-temperature limit of 
the retarded propagator is also causal (see Section 4.2.4 and Appendix C), this provides 
strong evidence that F{u,p,u) is causal in general. 



4.2.4 Zero-temperature limit 

As we move farther away from the horizon, we expect that near the boundary, the retarded 
propagator F approaches the AdS retarded propagator (which is expressed in terms of a 
Hankel function) for —p^\ sufficiently large. To simplify our discussion of the small u 

^This can be done in the WKB approximation. The proportionahty coefficient is fixed by solving 
F{u) for small u, such that it is normalized to 1 at u = 0, and constructing the interpolating function: 
i7r(cj + ie)'^u^/Tiu)/ ^/2^H^\{uj + ie)T[u)). 
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Re[F(a;,p,u)] Re[F(aj,p,u)] 




Figure 4: The exact F{u,p,u) (blue, solid) and its leading-order behavior (4.42) in the 
large u limit (red, dashed) with u = 0.5 and \p\ = 1. The plot on the right is a zoom at 

small values of u. 

limit we set p = 0. In this limit, and switching to Euclidean signature, the equation of 
motion for the massless scalar is 

<P'e - U'e - —<PE = 0. (4.45) 
u u 

The general solution to this equation is 

(t)E = CiuK2{2y/uUJE) + C2Ul2{2y/uUE)- (4.46) 

This is, of course, the general solution to the scalar wave equation in pure AdS, i.e. the 
zero-temperature limit of AdS-S. The dependence in (4.46) actually drops out, since 
'^E\fu = Ez/2. 

By allowing the range of z = ^/u/{'^^TH) to extend from to infinity, keeping the 
Bessel function which is well behaved in the interior, 

Fe{u,uje) = 2uuIK2{2^ue), u < 1, (4.47) 

and analytically continuing this result to Minkowski space [ue —i{uj + ie)), we obtain 

F{u, u) = ITIUUJ 2i/^^^(2v^a;), M<1, (4.48) 

where ifj^"* is a Hankel function of the first kind. In Figure 5 we plotted both the retarded 
propagator F and the zero-temperature propagator at u = 0.1, with the two curves almost 
perfectly overlapping (the two propagators differ for small values of uj). For comparison, 
we also plotted the finite and zero-temperature propagator for u = 0.5. 
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60 



-1 
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Figure 5: The zero (in red, dashed) and finite-temperature (in blue, sohd) 
bulk-to-boundary propagators at u=0.1 (left) and at u=0.5 (right), with p = 0. 

4.3 Thermal Feynman propagator and AdS-S circling rules 

Gibbons and Perry [29] noticed that the Hartle-Hawking definition of the Feynman prop- 
agator in a black hole background corresponds to the "Kruskal" vacuum. In other words 
it is periodic in imaginary time, and therefore can be identified with a thermal Green's 
function. 

To better understand this, we turn to the Rindler background example. Starting 
from 2-dimensional Minkowski space, ds"^ = —dt^ + dx'^ with the identifications x = 
exp(p) cosh(r),t = exp(p) sinh(r) one finds the 2-dimensional Rindler space metric ds"^ = 
exp(2p) [dp"^ — dr"^). The Feynman propagator corresponding to the Rindler vacuum is 
proportional to ln((Ap + Ar)(Ap — At)). On the other hand, the Feynman propagator 
corresponding to the Minkowski vacuum is proportional to ln((Ax+At)(Aa; — At)). When 
writing the Minkowski vacuum Feynman propagator in Rindler coordinates, one discovers 
that it is periodic in imaginary Rindler time: ln(cosh(r) — cosh(p)) + /(p). We arrive at 
the same answer when starting from the Rindler space retarded propagator in momentum 
space, in conjunction with the relation 



where /3 = 1/T = 2it. Lastly, we recall that the 2-dimensional Feynman propagator 
computed at finite temperature in imaginary-time formalism is ln(cos(r£;)— cosh(p)). Thus 
the Feynman propagator derived from (4.49) is a thermal Green's function by construction. 

Based on these observations, we proceed to define the Feynman Green's function 
corresponding to the "Kruskal" vacuum of the AdS-S black hole by 



Gf{E, P) = ReGniE, P) + i coth{f3E/2)lmGR{E, P), 



(4.49) 



Qf{uj,P,u) = ReQR{u,p,u) + i coth.{u7i)lmQ ji{u , p, u) 



(4.50) 
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where we used the fact that u = E / {2ttTh) ■ Since the retarded Green's function is analytic 
in the upper half-plane, (4.50) ensures that Qp will be a thermal Green's function, periodic 
in imaginary time. At this stage we can compute the other Green's functions by the usual 
relations: 

G' =Qf- Gr, G^ = Gf- Ga- (4.51) 

We can represent these bulk-to-boundary propagators using the same circling rules as 
before: the bulk-to-boundary propagator between circled and uncircled vertices is V~, 
between uncircled and circled vertices is and the propagator between two uncircled 
(or two circled) vertices is I^i;' (or Vp, respectively). 



5 Real-time finite-temperature AdS / CFT correlators 
5.1 The scalar two-point function 

In this section we review briefly the scalar 2-point function computation in real time 
and at finite temperature, following Son and Herzog [8]. We have already discussed 
the difference between Euclidean and Minkowski signature AdS-S geometries. We are 
interested in clarifying the issue of what bulk region must the tree-level gravity diagrams 
be integrated over. 



5.1.1 Minkowski Signature: the Schwinger-Keldysh Propagator 

From the work of Son and Herzog [8] we know that in real time, the gravity fields need to 
be specified on both time-like boundaries of the Penrose diagram. Therefore, one expects 
that one needs to integrate over the position of the bulk vertices over at least the R and L 
quadrants of the Penrose diagram. To obtain the result quoted in [8] , where the 2-point 
Schwinger-Keldysh field theory propagator is given as a boundary term, with the R and 
L contributions subtracted from each other, we will work with the bulk action 

S = N j y^(^^''d^(i)d^(t) + m^4>^ + interactions^ 

- iV y ^f^(^'"'d^(l)dy(t) + m^(l)^ + interactions^, (5.1) 

where iV is a supergravity normalization factor which includes the volume factor of S^, 
which was integrated over implicitly (to be specific, N = —N^/{167r^R^), where we recall 
that Nc is the number of colors in the dual A/'=4 super Yang-Mills theory, and R is the 
radius of S^). 
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In (5.1) the bulk right and bulk left contributions, from the respective boundaries to 
the event horizon, come with opposite signs ab initio. A similar action was used by Frolov 
and Martinez [11] in their construction of the action and Hamiltonian of eternal black 
holes in what they called "tilted foliation" . The Hamiltonian received two opposite sign 
contributions from the two causally disconnected regions, namely the R and L quadrants. 
Frolov and Martinez [11] pointed out that the Fock states in the R and L regions are akin 
to the doubling of the degrees of freedom (physical / doubler) encountered in the real-time 
formalism of finite-temperature field theories. Moreover, the structure of the Hamiltonian, 
a.s H = Hr — Hl, is also similar to the real-time Hamiltonian, which includes a physical 
particle Fock space contribution and a doubler Fock space contribution. 

The Kruskal coordinates tx and xk given in (4.5) are useful when considering global 
aspects of the spacetime and for defining notions such as incoming/outgoing and pos- 
itive/negative frequency modes. However, these coordinates are cumbersome when it 
comes to solving the equations of motion for and computing its action. For these tasks, 
it is simpler to return to the coordinate u, in terms of which the AdS boundaries are given 
by M = 0. However, u is only defined for one half of the extended AdS Schwarzschild space- 
time, so we need to define two such coordinates, ul and ur. ur is precisely our original 
u coordinate. We have 




It may seem that there should be a minus sign in the second equation since ul increases 
with increasing xk-, while ur decreases with increasing xk- However, Ir and ti flow 
in opposite directions with respect to tx- (tR flows roughly parallel to tx, while tz, is 
antiparallel.) 

Only the boundary terms of the quadratic action Sq play a role in the computation 
of the 2-point function. Since the boundary is perpendicular to the u direction, only the 
M-derivative terms, 

So = N j d'x j\uR^g--{d^(Py - N j d'^x j\uL^gg^-{du<Pf + ... (5.3) 
contribute to the boundary terms. Writing 

0(x,u) = j ^^e'P-^^{p,u), (5.4) 
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where = {E,P), and integrating by parts, we obtain the boundary terms^°: 

N f . 



(5.5) 



The 2-point functions are obtained by imposing the following boundary conditions on the 
scalar field [8]: 

i) (f>{p,u) is such that in the right quadrant, at the horizon, the positive energy modes 
are incoming and the negative energy modes are outgoing; the left quadrant modes are 
then uniquely determined by analyticity; 

ii) ) (j){p,u) approaches two distinct functions at the two boundaries: 

lim 4>{p,Ur) = 4>i{p), lim (f){p,UL) = hip)- (5.6) 
With this boundary prescription, the bulk scalar field in the R and L quadrants reads 

(f){p, Ua) = (t)h{p)Qha{,P, U), (5.7) 



where 



Q21 = 2^3- -lmF{-p,UR) 



' g2a;7r ^ 



^12 = -2^ g2^._i I^^(-P^^^) 



g2tJ7r 



/22 



F{-p, ul) - -^—-F{p, ul). (5.8) 



The indices a and b take the values 1 and 2, where 1 corresponds to the R quadrant, and 2 
the L quadrant. We have adopted the convention that, for example, Q12 is a propagator^^ 
which extends from a point on the R boundary to a point in the L bulk. By identifying 
F with the retarded bulk-to-boundary propagator, 

Fip,un)=gR, (5.9) 



^"^The horizon terms have been thrown away, based on the same considerations as in [7], namely the 
zero-temperature Hmit is correctly reproduced by the boundary terms at u^.i? = 0. 

^^Thc relationship between fk{u) in [8] and our F{u!,p,u) is complex conjugation. 

^^Strictly speaking, Qab are scalar bulk-to-boundary 2-point Green's functions, as indicated by the use 
of the letter Q. However, we will keep referring to Gab as scalar propagators, and we hope that this is not 
a source of confusion. 
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then the first equation in (5.8) is consistent with our definition of the Feynman bulk-to- 
boundary propagator (4.50) 

Gii = gF. (5.10) 

Plugging these solutions into the above action and varying with respect to 0i and 02 
yields a 2 x 2 matrix of 2-point functions Gab- It is important to recall that even though 
the scalar field and its action were real to begin with, the boundary conditions one had to 
impose at the horizon break reality. The on-shell action evaluated on the solution (5.7) 
is therefore complex and yields complex 2-point functions 

Gab{puP2) = -(-1)°+" , . (5.11) 

More concretely, the 2-point functions are 

(e^^" - l)F{p)d^F{p) + {e-^^-^ - l)F{-p)duF{-p) 



On = -2Ny^g 



uu ' 



(g2a;7r _ l)(e-2^^ - 1) 

Gi2 = -2iVv^^- ^^;^ _ ^ lFi-p)duFi-p) - Fip)duF{p) 



u=0 



it=0 



21 — —'-^127 

G22 = -Gil, (5.12) 

where we stripped away a momentum delta-function. This differs slightly from the 
result quoted in Son and Herzog [8] in that their retarded 2-point function is Or = 
—2Ny/gg'^^F{—p,u)duF{p,u)\u=o, whereas we find that 

Gnip) = -2N^g^^F{p,u)duF{p,u)\u=o. (5.13) 

However, since at the boundary F{p, u) is normalized to 1 and the term linear in u in 
F{p,u) is insensitive to the replacement of {uj,p) by {—u,—p), this difference is mostly 
cosmetic. Substituting (4.27), and using the small u expansions in (4.20), the retarded 
2-point function evaluates to^'^ 



Gr{p) = -AN{nTH)^R^ (^4 + C2 + 2A + 2c'2 + ^-^u + 2uci{l + + + divergent^ 
= -^-^ {cl + C2 + 2A + 24 + ^ ^ + 2wci (1 + z) + iu^ + divergent j , 



(5.14) 



^■^Wc have followed Son and Starinets by normalizing F{u = 0) = 1 and regularized Guip) by computing 
it at u = ub <C 1. Alternatively, we could have normalized F{u = ub) = 1, and found that equation 
(5.14) changes by the addition of — ^ ^2'^" {^^ — N^)^- 
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with the divergent terms given by 



divergent = ( ^ ^ cu + cA — + 2c2 \n{uB) 
\ Z J Ub 

= ^^^^-iu'-WYHuB). (5.15) 
Ub 

The coefficient A was computed in the hmit of small frequency in (4.29). The divergent 
terms can be dealt with either by subtracting the zero-temperature result for the retarded 
2-point correlator, or by employing boundary renormalization. As we will see in the next 
section, the finite-temperature CFT retarded 2-point can be obtained by analytically 
continuing the Euclidean result, yielding a natural expression in terms of a retarded and 
advanced supergravity mode propagator (a scalar mode in the case under consideration), 
just as we have obtained it from the R-L bulk action prescription, upon using energy- 
momentum conservation. 



5.1.2 Analytic continuation 

In the context of AdS/CFT, it has been pointed out by Gubser et al. [16] and by Iqbal 
and Liu [17, 18] that retarded CFT correlators can be obtained from Euclidean correlators 
by performing an analytic continuation. As before, only a boundary term in the classical 
action contributes to 2-point functions: 

^o = -N [ -^^^g^^FE{-k,u)d^FE{k,u)M-k)Mk) ■ (5.16) 

Since the integrand is singular as n — )■ 0, we placed the boundary at n = ^ 1. 

Varying the Euclidean action with respect to (p{—k) and (f){k) yields the 2-point func- 
tion: 

GE{k) = -N^g^'^FEi-k, u)duFE{k, u) . (5.17) 

u=0 

The analytic continuation from Euclidean to Minkowski space in (5.17) will produce either 
the retarded or the advanced propagators in the right-hand-side of (5.17), depending on 
the sign of ue- In particular, for ue > 0, this reduces to the previous expression found 
via the bulk action (5.1), namely N y/gg'^^F{p,u)duF{p,u)\u=o. 

5.2 Real-time finite-temperature 3-point functions 

Finally, we come to the main result of this paper, which is the prescription for computing 
CFT 3-point functions at finite temperature and in real time via supergravity diagrams. 
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Once we have the bulk-to-boundary propagators Qab, it is straightforward to compute 
3-point functions. These are obtained at tree level by computing the 0^ term (with 
coupling A say) in the action (5.1), plugging in the same solutions to the equation of 
motion for as in (5.7), and varying with respect to the boundary fields 0i and 02 as in 
(3.5): 

Gabc = -{-ir+'+'AN(^j^ dUR^galGblQcl - ^ dUL^Qa2Gb2Qc^. (5.18) 

where we recall that Qab were defined in (5.8). 

We would like to check if this expression for Gabc which we have obtained from the R 
minus L quadrant prescription for the AdS-S supergravity action (5.1) is correct. A first 
check would be to verify whether (5.18) obeys known identities for Schwinger-Keldysh 
3-point functions. In particular the KMS identities for the 3-point functions read 

G*abc = Gabc, (5.19) 

where 1 = 2 and 2 = 1, and we have assumed that a = (5/2. Using that the bulk-to- 
boundary propagators Qab have the property 

Qab = Qli, (5.20) 

it can be shown that, starting from (5.18), the identities (5.19) are indeed satisfied. 

Yet another test which (5.18) passes is to verify the largest time equation (3.15). On 
the gravity side we can think of this identity as arising from summing tree-level scalar 
field diagrams, with three vertices on the boundary of the R quadrant and one vertex in 
the bulk, with each vertex being either circled or uncircled for a total of 2^ diagrams. As 
expected in a causal theory, the sum is zero. 

Next, following [20] we define the "causal" (retarded) 3-point correlator to be given by 
(3.12). As shown by Kobes [20], the causal n-point functions are special, because they 
make contact with the n-point function computed in imaginary-time formalism. The 
causal real-time Green's functions are obtained from the imaginary-time Green's function 
by analytic continuation. 

When substituting the various Gabc in terms of bulk integrals in (3.12), given the 
multitude of terms that are added, we find that the final expression is surprisingly simple^^: 

GR{q,p;r) = AN6\p + q + r) [ duy/^F*{q)F*{p)F{r), (5.21) 

Jo 

similar conclusion was reached in [30] regarding the retarded 3-point function. 
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Figure 6: The retarded 3-point function from circling rules 



where F{r) is the retarded bulk-to-boundary propagator Qji. In hindsight, the retarded 
3-point correlator could only take the form in (5.21). That is because, by the circling 
rules, the retarded 3-point function can be written as a sum of R-quadrant tree-level 
gravity diagrams, with all vertices being either circled or uncircled with the exception of 
the vertex which has the largest time and which is uncircled. These are the diagrams 
shown in Figure 6. Their sum reduces to a single tree- level diagram, with two advanced 
and one retarded bulk-to-boundary propagator. This is the only possible answer for a 
causal theory when computing a tree-level retarded 3-point function. 

As another check on (5.21), we can take the zero-temperature limit. In this limit, the 
retarded/advanced bulk-to-boundary propagators of a massless scalar field in the AdS-S 
background approach the AdS retarded/advanced propagators (given in Appendix C), as 
discussed in Section 4.2.3. The integration limits extend in the zero-temperature limit to 
(0, oo), and the integral in (5.21) reproduces (2.22). We have also verified that (5.21) can 
be obtained by analytically continuing the Euclidean space 3-point function. 



A Momentum-space 2-point CFT correlators at T=0 

We now want to find the momentum-space expression of all 2 and 3-point functions given 
above. The 2-point correlator in momentum space, Euclidean signature is 

GEik) = -. '^''"f.."t"'> -)^-'^-, (A.1) 
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where k"^ = E"^ + k"^, and the Fourier-transform integral was regularized by dimensional 
regularization. Note that if the conformal dimension A is a positive integer greater or 
equal to 2, then the 2-point will contain logs. Not surprisingly, this result coincides with 
the AdS/CFT computation. 

In Minkowski signature, the story is a bit more complex. The Fourier transform of 
the scalar field 2-point Green's functions yields: 

G± = -2Tn9{±E)6{E^ - P^), (A.2) 

Gr{p) = - g, , Ga{p) = GUp), (A.3) 

-{E + 2e)2 + P2 

Gf{p) = ^ . , (A.4) 

-E^ + -te 

while the momentum-space expressions of the non-time ordered, retarded, advanced and 
time-ordered 2-point scalar field operator correlators are 

= -^^^^^A±E-\P\)iE^-PY^^ (A.5) 

Gr;A>2{p) = J^^^^^,{E'-PY-'H-{E + ^er + P')e') (A.6) 



(^2 _ p2)A-2 i^^d _g2^ p2|^2) _ ^^Q^g2 _ p^^^g^^E) 



(2^-ir(A)) 

Ga{p) = GUp) (A.7) 

~ ■'■) ^ Z72 d2\A-2/ Tp2 , d2|,2^ 



GF;A>2ip) = (2A^ir(A))2 (-^ - PT [H\ -E' + PY) - ^^e{E' - p')) , (a.s) 

where A = 1/e was used as a UV cut-off to regularize some of the integrals, and A was 
assumed to be an integer, for concreteness. 

Note that the analytic continuation of the Euclidean signature 2-point function, with 
E — > —i{E ± ie) yields the retarded/advanced 2-point correlators. 



B Momentum-space real-time 3-point correlators at 

r=o 

Here we present the intermediate steps leading to the results quoted in Section 2.2. We 
begin with the Euclidean 3-point correlator in position space, and Fourier transform. 
Since the integral is not convergent in c? = 4, we use dimensional regularization, with 
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d = 4:-2e: 



*^12'^23'^31 



^ i=i ^ 



,-1 Jo 



(SiS2 + S2S3 + S^Si)i 



4 



• , >/vj / V Uin2U3(— + — + — ' 



Vi=i^o 



r(35- 



M— - 2'' ^ ' 



(B.l) 



which is the expression we gave in (2.16). 

Next, we compute the Fourier transform of (2.21) first in the kinematics 1,2 incoming 
and 3 outgoing, that is i?i < 0, i?2 < 0, -Ea > 0. The steps taken next follow closely the 
path used for the Euclidean 3-point correlator. We could have chosen to regularize the 
integrals using dimensional regularization, as we did before. However, here we chose a 
different regularization to show how the momentum-space correlator can be obtained via 
a 5-dimensional AdS integral. At each step in our manipulations we will make sure that 
we appropriately insert convergence factors, such that the integrals are well defined. 

To begin, we use again Schwinger parameters to write the denominators of the first 
term in (2.21), which is the only term that contributes in the chosen kinematics. For 
example. 

The integrals over the position-space coordinates are of the type 

dze'"'^ = /^e^^§''(")"/^ a e R. (B.3) 
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The intermediate result, where we still have to perform the integral over the three 
Schwinger parameters, is very similar to what we have encountered in Euclidean sig- 
nature: 



exp (-^(uiPi + UaPa + ^sPs) ) ' (B-4) 



i-l 

1 / i 



_|_ J_ _|_ J_)35-2 



where p\ = —Ef + etc. The integrand in (B.4) is properly regularized for both Ui = 
and for large values of ut. To deal with the + — + — ) factor, we introduce another 

Schwinger-type parameter, but in a slightly different fashion than previously 



1 2e-^ 



(a + ie')'" r(n 
The three Schwinger parameter Ui integrals are of the type 



POO 

/ dzz2"-ie^("+^^"')^', aeR. (B.5) 
Jo 



(B.6) 

In what follows we set the convergence factor e = 0. Thus we find that the retarded 
3-point function in momentum space can be written as 

(27r)8 ^ dz 

Gn{PuP2;Ps) = 23A-ar(A/2)3r(3A/2-2) ^^^^+^^+^^Vo 



z\-El + - 2e) ^-'K^^2{z\ -El +f^-ie 



z'i-El + pi + ie)^~^K^^2{zJ-El +f^ + ie 



(B.7) 



A convergence factor e"*^ ^ was also dropped since the integrand is convergent at large z. 
To simplify notation we removed the tilde from e. 

From the largest time equation we learn that the causal n-point function is real: 

Gr{{^}) = GU{^}) (B.8) 
which implies for the momentum-space causal retarded n-point function 

Gni{p}) = GrU-p}). (B.9) 
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In particular, this gives us the causal retarded 3-point function for the energy configuration 
-El 2 > and < 0, and leads us to the following expression for Gji{pi,p2;P3) for general 
kinematics 

Gn{v..V2;v.) = 23A-6r(A/2)3r(3A/2-2) ^ (P1+P2+P3) | 

z\-{E, - ief+f,)^-'K^.2{z^-{E2 - lef+f^) 
z'{-{E, + ie)f + f,)^~'K^.,{z^-{E, + leY + pf) . 

(B.IO) 



C The retarded bulk-to-boundary scalar propagator 
in AdS space in Poincare coordinates 

In this section we briefly discuss the causal properties of the retarded/advanced bulk-to- 
boundary scalar propagator in 5-dimensional AdS 



Gr = C z\-{E + ie)f + f)^-^K^^,{z^-{E + ieY+f), (C.l) 

where C is a real- valued normalization constant. The simpler route to constructing the 
retarded propagator in position space is to first compute the Feynman propagator 

gF{EJ,u)=gR{E,p,u)ii E>Q, and m) = m) if E < (C.2) 

which leads to 



Gf = Cz\-E^+f- ie)"^K^_2{z^-E^+p'-ie). (C.3) 

Then the Fourier transform of the Feynman propagator can be computed with the help 
of (B.6) 

This is the expression we might have gotten starting from the Feynman bulk-to-boundary 
scalar propagator computed in Euclidean AdS, after performing the usual Wick rotation 
and analytic continuation. 
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Hence the retarded propagator in position space obtained from 

QR = o{t){gp + g*p) (C.5) 

is causal, and has support only on the forward hght cone. 

D The retarded 3-point momentum-space correlator 
with A=2 and at r=0 

In this appendix we give a closed-form analytic answer for the retarded 3-point CFT 
correlator, of three scalar operators with conformal dimension A = 2. Again we begin 
from the position-space (2.14), and we perform the Fourier transform as follows: 

i^GR{pi,p2;P3) = n / d^Xie-'P'''H^GR{xi,X2,X3) 

■ /'.---■■7'™--/'--/-/S=- - - 



dzi dz2 



2Txi[zi — ie) J _^ 2ni{z2 — «e) 



,X2;X3) 

'4^ / ^4^ / j4:„ 



d Xi2 d X23 / d X31 



gi(£;i+^l+n)ti2g-i(Pl+-ft:)xi2g-i{£;3-n)t23g»(J'3--f^)^23g*{f^+^2)t31g-«^^^31 

1 1 



c.c I I — — c.c. 



xi2 + ieti2 - iet23 J + iehi 

_j_g«{-Bi -zi +^)ti2 g-«( A +K)xi2 g-i(-E3+22-r2)i23 g*(-P3 -K)x2-i ^int3i ^-iKx-ii 
1 1 \ ^ 1 



— C.C I I — — c.c. 



xf2 - "^£^12 2:31 + ietsi J \X23- iet23 



Next substitute the Fourier transform of the functions depending on X12, X23, X31, and use 
that these are of the type 

9i±E)6iE^ - P2) = ^^'^'J^^ (D.2) 



(D.l) 
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to perform the integrals over Zi, and f2. The remaining d^K integral can be simplified 
by going to a Lorentz frame where P3 is vanishing: 



{4:71 



3\3 



(27rz)2 
1 



1 



-2\K\ 



|Pi + K\ + \K\ + E^ + E^ + te {\K\ + E3 + lef - \K\' 
-2\K\ 



]Pi + K\-\K\ + E^ + E^ + ie 
1 



+ 



■\K\+E-i + ieY-\K\^ 
-2\K\ 



|Pi + K\ + \K\ -Ei+te i\K\ +E3 + leY - \K\^ 
1 -2\K\ 



-|Pi + K\ - \K\ -Ei + te {-\K\ +E3 + leY - \K\^ ) \K\^\Px + K\ E3 + % 



le 



(D.3) 



The 3-dimensional integral over K is performed by using spherical coordinates, with Pi 
aligned with the z-axis. The integral over is trivial and yields a factor of 2ti. The 
remaining integrals over \K\ and 9 are re-expressed as integrals over \K\ and \Pi + K\. 
The Jacobian of this change of variable is 



(9(|i?|,cos^) 



d{\Kl\P^ + K\) 



\Pi + K\ 
\Pi\\K\ 



(D.4) 



where |i^'-|-Pi| is integrated from Hii'l — |Pi|| to |i^| + |Pi|. Lastly the integral over |i^+Pi| 
is easily evaluated to a log, and the final expression of the retarded 3-point function is 



i^GR{pi,p2;p3) 



E3 + ie 
1 



d\K\ X 
1 



In 



2\K\ + \Pi\ + El + E3 + ie 



P3 + 2\K\ + te I \K\ - |Pi| I + \K\ +Ei + E, + ie 



In 



-2\K\ - \Pi\ + El + E3 + ie 
E3-2\K\+ie -\\K\ - lAll - + ^i + ^3 + ^e 

2\K\ + \Pi\-Ei + ie 



E3 



le 



In 



E^ + 2\K\+ie Ili^^l-lPill + li^l-Pi + ze 
-2\K\-\Pi\-Ei + ie 



In 



E3 - 2\K\ + le -\\K\ - \Pi\\ - \K\ -Ei + tej } 



(D.5) 
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With some effort, the integral over can be performed analytically: 



(47r)3(27r) 



{2my\Pi\{Es + ie) 



2Li2 



-Eo 



\Pi\+ie 



Ei + \Pi\ 



2Lu 



E2 + \Pi\-ie 
-El + \Pi\+ie 



-In' 
2 



-E2 + \Pi\+ie\ 1 , 2 

In 

2 



E, 



El + \Pi 
+ \Pi\ 



^2 + |Pi| 



le 



El + \Pi 



le 



-El + I Pi 
-E2 + |Pi| +ie 
-El + iPil + «e 



(D.6) 



where we recall that we went to a special Lorentz frame such that P3 = 0. We have 
checked that in this frame, our result coincides with eqn. (28) of [31]. The final form of 
the Fourier-transformed retarded 3-point function, with ^3 being the momentum flowing 
in the vertex with the largest time, is obtained by covariantizing the previous expression 
with the help of the following relations 



l^il 



El 



{pmY_ 



El 



{P2Pzf 

pI 



(D.7) 



E The Heun's Function 



Heun functions arise as a generalization of the hypergeometric function, and are defined 
as solutions of the Fuchsian differential equation [32, 33, 34]: 



7 



5 



dH 



af3z — q 



z z — 1 z — d J dz z{z — l){z — d) 



H = 0, aP,-f,6,eeC. (E.l) 



This differential equation has 4 regular singular points: 0, 1, li, 00. Regularity at infinity 
is insured provided that 

e = a + {3--f-5 + l. (E.2) 

The characteristic exponents at the singular points are: 0, 1 — 7 for z = 0; 0, 1 — 5 for 
z = 1; 0,1 — e for z = d, and a, /3 for 2; = 00. 

The Heun's function is further defined to be normalized to 1 at z=0: 



Hl{d,q,a,P,-f,6] z) 



1: 



z=0 



d 

dz 



Hl{d,q,a,P,-f,6]z] 



2=0 



'yd 



(E.3) 



In the vicinity oi z = 0, the Heun's function is analytic and can be expressed in terms 
of a locally convergent series expansion X]fc>oCfc^'^, where the coefficients Ck are obtained 
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recursively 



{k + l){k + j)dCk+l 



ki {k + J + 5 - l)d + {k + + e - 1) ) + q 



Ck 



+ {k + a- l){k 
c_i = 0,co = 1. 



(3 - l)ck-i = 



(E.4) 



The series expansion breaks down when k is a negative integer, and the Heun's function 
will exhibit logarithmic behaviour. This case will be further investigated, as it is relevant 
for the physical problem we are considering. 

Based on the characteristic exponents at 2; = 0, there are two independent solutions of 
the Heun differential equation. One is Hl{d, q, a, (3, 7, 6; z), and the other is z^~"'Hl{d, q — 
(7 — l){6d + e), /3 — 7 + 1, a — 7 + 1, 2 — 7, 5; z). Notice that for the latter function, 
having 7 be a negative integer poses no special problems. Then, we can search for a 
particular combination of these two independent solutions such that the poles in 7 cancel. 
In particular, for 7 = —1, the linear combination 



lim Hl(d, q, a, f3, 7, S, z 
7— 1 ' 



C2 



7+1 

where C2 is obtained from (E.4) 



z^~'<Hl{d, g-(7-l)(M+e), /3-7+1, a-7+l, 2-7, ^; z) 

(E.5) 



(7 + + 7 + e + g ) ^ - a/3 



(E.6) 



is well defined near z = Q. In (E.5), the limit is taken with c?, a, /3, 5, 2; held fixed and with 
e treated as a function of 7 accordining to its defining expression e(7) = a + — 7 — 5+1. 
The leading term from the second Heun function will cancel the 2:^/(7 + l)-terms from 
the series expansion of the first Heun's function. The logarithmic behaviour of (E.5) 
arises from the prefactor z^~'^ /{;~^ + 1), written as exp[(2 — (7 + l))ln(z)]/(7 + 1) ~ 
^2(1/(7 + 1) -ln(z) + 0(7 + 1)). 

The series 'Ylik>o^kZ^ converges for \z\ < min(l, \d\). However, through analytic con- 
tinuation it is possible to extend the Heun's function on the complex plane. There are 
192 known local solutions of the Heun differential equation [34] which allow the Heun 
function to be extended in different regions of parameter space. For example. 



Hl{d, q, a, /3, 7, 6; z) 



i7/(c/,g-(5-l)7rf,/3-(5 + l,a-5 + l,7,2-5;^). (E.7) 



Or, one can instead choose to express the solution of the Heun differential equation as an 
analytic expansion around z = 1 instead of z = 0: 



HI{1 — d,—q + af3, a, 13, 5, 7; 1 — z) 



(E.8) 
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is the solution with characteristic exponent 0, and normahzed to 1, while 

{z-iy-^Hl{l-d, -q+{6-l)-fd+{(3-6+l){a-6+l), 13-6+1, a-6+1, 2-6, 7; 1-z) (E.9) 

is the solution with charactersitic exponent 1 — 6. The complete list of the 192 different 
local solutions of the Heun function can be found in [34] . 
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